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Stability control in laser is still an emerging field of research. In this paper the dynamics of 
External cavity semiconductor lasers (ECSLs) is widely studied applying the methods of chaos 
physics. The stability is analyzed through plotting the Lyapunov exponent spectra, bifurcation 
diagrams and time series. The results of the study show that the rich nonlinear dynamics of the 
electric field intensity (| E p) includes saddle node bifurcations, periodic, quasi periodic, chaotic 
and regular pulse packages. The issue of finding the conditions for creating stable domains has been 
studied. By considering the dynamical pumping current system coupled with laser, a method for 
the creation of the stable domain has been introduced. 
Usage: 

PACS numbers: 42.65.Sf, 42.55.Px, 05.45.Xt,05.45. Jn, 05.45.Vx 



I. INTRODUCTION 

The research on the laser dynamics has attracted a lot 
of interest for the past several years. Current investiga- 
tions in this field include the study of different routes to 
chaos in lasers and the control of chaos in such systems. 
Most of these works are concentrated on semiconductor 
lasers with respect to their importance in optical com- 
munications. Chaotic behavior in semiconductor lasers 
has been investigated under various physical conditions 
such as external optical injection, and external optical 
feedback 

Our study is motivated from different points of view, in 
this study dimensionless ECSLs are investigated [3] ■ This 
class of dynamical systems can successfully be modeled 
by delay differential equations. Where the study has pro- 
vided a realistic background for the investigation of time 
delay systems, which is very difficult as a problem in the 
infinite-dimensional space [5|. This allows us to present 
a consistent overall picture of the dynamics at the same 
time. In this study parameters such as feedback phase, 
feedback strength, and pump current are used as control 
parameters to obtain different dynamical regimes. 
Also the key idea underlying most chaos controlling 
schemes such as ECSLs is to take advantage of the un- 
stable steady states and unstable periodic orbits of the 
system They are embedded in the chaotic attrac- 

tor characterizing the dynamics in phase space. Feedback 
schemes require real-time measurement of the state of the 
system and processing of a feedback signal 9, 1^. The 
location of the unstable fixed-point must be determined 



before control is initiated. The standard control method 
has the disadvantage of ha ying a small region of conver- 
gence around the fix point |lll - ll3j . 

The aim of the present paper is studying the ECSLs dy- 
namics by the Lyapunov exponent diagram as the ana- 
lyzing tool and proposing a control method for the Laser. 
Regarding the pump current as one of the control param- 
eter of the system, a discrete time dynamical system is 
built and coupled with Laser. 

The rest of the paper has been organized as follows: Sec- 
tion 2 describes Lang and Kobayashi model for ECSLs. 
The method of the control and analyzing the results are 
proposed in Section 3 and 4. In section 5, the stable do- 
main of the introduced model before and after applying 
the control method is studied via Lypaunov exponents 
and bifurcation diagrams. Section 6 concludes the let- 
ter. The paper ends with an Appendix which contains 
algebraic calculations of the controller. 



II. THE LASER MODEL 

In the early 1980s, Lang-Kobayashi (LK) proposed 
model semiconductor lasers. The LK equations are model 
equations that have been used extensively in the past to 
describe a semiconductor laser subject to feedback from 
an external cavity [3 [l^. LK equations are a sort 
of delay differential equations (DDEs) with an infinite- 
dimensional phase space |3|]. For the (complex) electric 
field E and inversion N, we write the LK equations as the 
dimensionless and compact set of equations [3] : 
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Parameters describe the line width enhancement factor 
a, the feedback strength rj, the 27r-periodic feedback 
phase Cp, the ratio between carrier and photon hfetime T 
and the pump current P. In these equations, the time is 
normahzed to the cavity photon hfetime (Ips) [l6|. The 
external round trip time r is also normalized to the pho- 
ton lifetime. The remaining parameters are held fixed at 
T = 1710, P = 0.8, r = 70, a = 5.0 [l^. In this paper, 
the stability of an electric field intensity | ^ is stud- 
ied versus the feedback phase change Cp and feedback 
strength change r]. 



III. THE CONTROL METHOD 

Most recently, a different method of control has been 
shown to be successful in the experimental control of 
chaos. Based on: 

• Determination of the stable and unstable directions 
in the Poincare section. 

• Self-controlling feedback procedure. 

• Introduction of small modulation of a control pa- 
rameter. 

• Knowledge of a prescribed goal dynamics. 

The first two methods are usually called feedback meth- 
ods, while the 3rd and last are called non-feedback meth- 
ods [3| • Although control of chaos by small modu- 
lations has not been proved in general [20[, this method 
involves the on chaotic behavior generation of an error 
signal from the difference between the output signal and 
its value at an earlier time. A great virtue of this method 
is that it does not require knowledge of more than one 
variable. The method should be very useful for apply- 
ing in fast systems to control a chaotic nonlinear circuit. 
Until now, different parameters such as modulating cur- 
rent [111, , modulating voltage [l^ , and pump power 
variation response to a continuous error signal genera- 
tor have been used to modify the laser output. In 
this study, "pump current" has been selected as a tool to 
be implemented in the control method. By considering 
the Laser and the control system as a two dimensional 
dynamical system, the simple model for controlling the 
stability of Laser is introduced as follows [Isl] : 
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The setup is shown schematically in Fig. 1. The most 
significant parts of the setup include the laser pump and 



the method whereby its output is modulated. The pump 
power is controlled by the use of an electric circuit. Fig- 
ure 2 shows the desired controller. A photodiode is used 
to convert a sample of laser light to electric voltage "K,". 
This voltage is connected to non-inverting input of the 
comparators "C/i" and "C/2". In the stable condition, 
Vs is less than Vi and more than V2. Thus, the output 
voltage of Ui and U2 is in the low and high conditions, 
respectively. Consequently, Qi and Q2 are off and Q3 
is on. As a result the relays 1 and 2 are off and is 
connected to the input of the Pm+i function generator. 
When the unstable condition occurs and the laser inten- 
sity is increased, the voltage Vs becomes more than Vi 
and then the output voltage of Ui changes to high. Due 
to high condition in the output of Ui, the transistor Qi 
and Q2 conduct. At this stage the relay 2 disconnects 
P„j from the input of the Pm+i function generator and 
simultaneously relay 1 connects the Pm+i function gen- 
erator output to its input. The feedback from the out- 
put to the Pm+i function generator input continuously 
decreases the Pm+i, and consequently decreases photo- 
diode output voltage. Finally, when the output voltage 
of the photodiode becomes less than Vi, Qi and Q2 be- 
come off. In this condition relay 1 disconnects but relay 
2 is still on. And the capacitor "C" at the input of Pm+i 
function generator holds the last Pm+i- After a while, 
if laser intensity decreases more, and "Fs" becomes less 
than V2, then, the output voltage of U2 changes to zero 
and then changes to off. As a result, relay 2 becomes 
off and connects the Pm to the input of the Pm+i func- 
tion generator. Finally, Pm+i increases and the system 
changes to normal condition. 

In this study, arbitrary chaotic regions of the system dy- 
namical behavior, have been chosen to undergo the pro- 
posed control method [20l - [2^ . 



IV. STABILITY ANALYSIS 
A. Lyapunov exponent spectrum 

Lyapunov exponents and entropy measures, can be 
considered as "dynamic" measures of attractors complex- 
ity and are called "time average" [2^. The Lyapunov 
exponent A is useful for distinguishing various orbits. 
Lyapunov Exponents quantify sensitivity of the system 
to initial conditions and give a measure of predictabil- 
ity. The Lyapunov exponents are a measure of the rate 
at which the trajectories separate one from another. A 
negative exponent implies that the orbits approach to a 
common fixed point. A zero exponent means that the or- 
bits maintain their relative positions; they are on a stable 
attractor. Finally, a positive exponent implies that the 
orbits are on a chaotic attractor, so the presence of a pos- 
itive Lyapunov exponent indicates chaos. The Lyapunov 
exponents are defined as follows: 

Consider two nearest neighboring points in phase space 
at time and t, with the distances of the points in the ith 
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direction ||fai(0)|| and ||5a;i(<)||, respectively. The Lya- 
punov exponent is then defined by the average growth 
rate A.; of the initial distance, 



A,^lin.ilnf^. 

t^oo t \\SXi{0)\\ 



(4) 



The existence of a positive lyapunov exponent is the indi- 
cator of chaos showing neighboring points with infinitesi- 
mal differences at the initial state abruptly separate from 
each other in the ith direction \2f\. Using the algorithm 
of Wolf , the lyapunov exponent is calculated versus 
a given control parameter. Then, the value of the control 
parameter increases a little and the Lyapunov exponent 
is calculated for the new control parameter. By contin- 
uing this procedure Lyapunov exponent spectrum of the 
system is plotted versus the control parameter. 



B. Bifurcation diagrams 

Bifurcation means a qualitative change in the dynami- 
cal behavior of a system when a parameter of the system 
is varied [26j . A bifurcation diagram provides a useful in- 
sight into the transition between different types of motion 
that can occur as one parameter of the system alters. It 
enables one to study the behavior of the system on a wide 
range of an interested control parameter. In this paper 
the dynamical behavior of the system is studied through 
plotting the bifurcation diagrams of the of intensity | | ^ 
versus feedback phase Cp and feedback strength change 
77, as control parameters. This procedure continued by 
increasing the control parameter and the new resulting 
points were plotted in the bifurcation diagram versus the 
new control parameter. 



V. RESULTS AND DISCUSSIONS 

A. Introducing the dynamics of the master laser 

The dynamical behavior of the master laser as a func- 
tion of the control parameter, Cp, can be divided into 
two areas: 

• low feedback strengths 

• high feedback strengths 

In the low feedback strengths, such as; r]m= 0.0455 [l6| . 
the output of laser contains chaotic (CH) [2^, periodic 
(P1,P2,..) and quasi periodic (QP) HO], behaviors, re- 
spectively. Fig. 3(a) shows that the dynamical behavior 
of the output intensity for master laser is CH, P7, PI and 
QP. In the high feedback strengths, such as; rjm— 0.135, 
the laser output contains regular pulse package (RPP), 
PI, and QP behaviors which is confirmed in Fig. 

4. As it is shown, the dynamical behavior of the output 
intensity for master laser is RPP, QP and PI. For under- 
standing the general faces of the ECSLs one could refer 



to Fig. 5 , where the bifurcation curve of the laser output 
in terms of feedback strength 77, in a constant feedback 
phase Cp is presented. As it is shown, the bifurcation 
curve contains periodic, QP, CH, and RPP behaviors. 
The maximum Lyapunov exponents are presented in Fig. 
3(b) to verify the corresponding characteristics. As can 
be observed, the positive and negative values in Lya- 
punov exponent spectrum verify the aperiodic and pe- 
riodic behaviors in bifurcation diagram. In addition, ex- 
tended numerical analysis on the ECSLs dynamics has 
been carried out on several control parameters. The 
system exhibits extremely nonlinear behaviors when the 
control parameters are varied. 



B. Applying the chaos control method 

In order to provide a fair justification on employing a 
periodic perturbation as the control technique, it is first 
necessary to consider two facts. First, the similar chaos 
control methods should be evaluated, and their practica- 
bility and associated advantages should be clarified. Sec- 
ond, we address this question that: can such a method 
be easily implemented in the real world applications? In 
fact, a method should be chosen so as to be best suited for 
an experimental implementation. In this way, as it will 
be discussed in the following, the pump current sources 
will be able to satisfy both highlighted remarks. 
Before evaluating the potential chaos control methods 
for an ECSLs system, we should consider the parameters 
that can be implemented in the control technique. It is 
clear that, to control the chaotic oscillations of a dynam- 
ical system, one of the control parameters of the system 
should be perturbed or modified in order to reach the 
regular desirable behavior. However, in the application 
phase , we are restricted to apply changes only to the ex- 
ternal forcing elements such as feedback phase, feedback 
strength that are determined by the type of the applica- 
tion, and the pump current which has been chosen in the 
present study. 

The results are depicted in Figs. 6-9. Under the collec- 
tive system parameters, every category corresponds to 
original system and controlled system; see Figs. 6(a), 
6(b) and Figs. 7(a), 7(b). Also the control method has 
been tested through Lyapunov exponent diagrams; see 
Figs. 6(c) and 7(c). The dashed line represents the orig- 
inal system, while the solid line represents the controlled 
system. This figure indicates a significant abatement of 
the Lyapunov exponents from positive values to negative 
ones indicating that stable dynamics can be achieved af- 
ter the proposed technique is engaged. The results have 
also been confirmed by plotting the intensity \ E versus 
time in a certain value of the feedback phase Cp before 
and after the control process. Respectively Figs. 8(a) 
and 9(a) show before control process and Figs. 8(b) and 
9(b) depict after control process. 
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VI. CONCLUSION 

Our main contribution in this paper is to develop the 
simple method for study and control the laser dynam- 
ics. In this paper, the concept of feedback chaos control 
has been re-defined by taking into account the control 
parameter as a variable in the time and is changed by 
using another chaotic map, for which a new and effective 
control scheme has been presented. The result of the 
present study broadens our understanding of the com- 
plex dynamics of laser, and also helps us to control its 
nonlinear dynamics. 



ACKNOWLEDGMENTS 

The author expresses his sincere thanks to Professor 
A. Jafari for his invaluable advice and encouragement. 



Appendix A: Detail of control circuit 

The circuits in the Figs. 10, 11 and 12 are used to 
generate Pm+i Function generator. Figure 11 shows the 
multiplier circuit. This circiiit multiplies (1 — _P,„) by 
(1 — Pfn) and makes the denominator of the Pm+i func- 
tion generator. The derived equations for Fig. 11 are as 
follows: 

Vb = -Vd = -VrLni^ + 1) 



(Al) 



The op amp U5 adds the two equal voltages Vb to its 
inverting input. Then, the voltage 14 at the output f/5 
is: 



Vc = VT{Ln{^ + l)+Ln{ ^" 



RIs 



1) 



= VrLn 



(i 



Va 



VRIs 



(A2) 



and the voltage at the output of Uq is: 

Vd = -RIdi = -RIs {e^^^ - 1) 



y2 
« - 2V 

RIs 



(A3) 



The op amp Uj adds the voltages Vd and Va , and there- 
fore, cancels "-2Va" at its output. Then, the voltage at 



Vg will be as follows: 



y2 

V — ° 
'-RIs' 



(A4) 

By considering R = IMfl and Is = 1/isec, the equation 
{A4:) can be written as: 



Vg = V^ = (1 - Pmf 



(A5) 



In Fig. 12, a part of multiplier circuit including op amp, 
U2, U4, U5 and Ua is used as a divider. In this figure, 
the voltages at points a and b are respectively APm and 
— (1 — Pto)^ , therefore, the voltage at point c is as below: 



V, = Vt{ Ln& + 1) - in(-L^ + 1) 



RIs 



VrLn 



RIs 



4Pm + RIs 
{l-Pl) + RIs 



(A6) 



Assuming APm » Rh and (1 - PmY > Rh, the equa- 
tion {A6) can be written as below: 



Vc = VrLn 

and then Vd is: 

Vd = -RIs{e^ - 1) 
Finally Vf is: 

Vf 



iP,,:. 



-RIs 



AP 



(1 - Pmy 



AP 



{1-Pmr' 



(A7) 



(A8) 



(A9) 
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